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ABSTRACT 

The realization that the total energy of GRBs is correlated with their jet break angles moti- 
vates the search for a similar relation between the peak luminosity, L, and the jet break angles, 
L oc 8~ 2 . Such a relation implies that the GRB luminosity function determines the angular 
distribution. We re-derive the GRB luminosity function using the BATSE peak flux distribution 
and compare the predicted distribution with the observed redshift distribution. The luminosity 
function can be approximated by a broken power law with a break peak luminosity of 4.4 x 10 51 
erg/sec, a typical jet angle of 0.12 rad and a local GRB rate of 0.44 hg 5 Gpc yr . The angular 
distribution implied by L oc 9~ 2 agrees well with the observed one, and implies a correction factor 
to the local rate due to beaming of 75 ± 25 (instead of 500 as commonly used). The inferred 
overall local GRB rate is 33 ± 11 /i| 5 Gpc~ 3 yr _1 . The luminosity function and angle distribution 
obtained within the universal structured jet model, where the angular distribution is essentially 
oc 9 and hence the luminosity function must be oc L~ 2 , deviate from the observations at low 
peak fluxes and, correspondingly, at large angles. The corresponding correction factor for the 
universal structure jet is ~ 20 ± 10. 

Subject headings: cosmology:observations-gamma rays:bursts 



1. Introduction 

In spite of great progress, one of the missing 
links in our understanding of GRBs is a knowledge 
of the luminosity function and the rate of GRBs. 
The realization that GRB are beamed raised the 
question of what is the angular distribution. The 
discovery of Energy - angle relation (Frail et al., 
2001; Panaitescu & Kumar 2001) suggested that 
the luminosity function and the angular distribu- 
tions should be also related. 

The number of GRBs with an observed after- 
glow is still rather limited. Redshift has been mea- 
sured only for a fraction of these bursts and the 
selection effects that arise in this sample are not 
clear. Therefore, at present we cannot derive di- 
rectly the GRB luminosity function. A possible 
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method to derive the luminosity function makes 
use of a luminosity indicator. Several such indica- 
tors have been suggested (Norris, Marani & Bon- 
nell 2000; Fenimore & Ramirez- Ruiz 2000) but the 
robustness of those indicators is not yet clear. In 
a recent paper Schmidt (2003) shows that a lu- 
minosity function obtained using luminosity indi- 
cators fits the observations very poorly. An al- 
ternative is to consider the observed peak flux 
distribution and fit possible luminosity functions 
and burst rate distributions (Piran 1992, Cohen 
& Piran 1995, Fenimore & Bloom 1995, Loredo & 
Wasserman 1995, Horack & Hakkila 1997, Loredo 
& Wasserman 1998, Piran 1999, Schmidt 1999, 
Schmidt 2001, Sethi & Bhargavi 2001) This can be 
done by obtaining a best fit for the whole distribu- 
tion or simply by just fitting the lowest moment, 
namely, (V/V m ax / ■ 

Evidence of jetted GRBs arises from long term 
radio observations (Waxman, Kulkarni & Frail 
1998) and from observations of achromatic breaks 
in the afterglow light curves (Rhoads 1997, Sari, 
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Piran & Halpern 1999). However the structure of 
these jets is still an open question. The two lead- 
ing models are (1) the uniform jet model, where 
the energy per solid angle, e, is roughly constant 
within some finite opening angle, 9, and sharply 
drops outside of 9, and (2) the universal structured 
jet (US J) model (Rossi, Lazzati & Rees 2002), 
where all GRB jets are intrinsically identical, and 
e drops as the inverse square of angle from the jet 
axis. Within the uniform jet model, the observed 
break corresponds to the jet opening angle. Frail 
et al. (2001) and Panaitescu and Kumar (2001) 
have estimated the opening angles 9 for several 
GRBs with known redshifts. They find that the 
total gamma-ray energy release, when corrected 
for beaming as inferred from the afterglow light 
curves, is clustered. A recent analysis (Bloom, 
Frail & Kulkarni, 2003) on a larger sample con- 
firmed this clustering around ~ 1.3 x 10 51 ergs. 
Within the USJ model the jet break corresponds 
to the viewing angle and the energy-angle relation 
given above implies that the luminosity distribu- 
tion within the jet is proportional to 9~ 2 (Rossi, 
Lazzati & Rees 2002) and therefore the luminos- 
ity function has the form $(£) oc L~ 2 (Perna, 
Sari & Frail 2003). Frail et al. (2001) also derive 
the observed 9 distribution. Taking into account 
the fact that for every observed burst there are 
= (9 2 /2)~ 1 that arc not observed they derived 
the true 9 distribution and used it to estimate the 
average "beaming factor", (f- 1 ) ~ 500. Then 
they multiply this factor times the local (isotropic 
estimated) GRB rate derived by Schmidt (2001) 
and obtain the true local rate of GRB. 

We aim here at obtaining the combined L-9 
distribution function, <&(L,6) in the uniform jet 
model. The present data is not sufficient to accu- 
rately constrain $(L,#). We therefore must make 
some assumptions. First we estimate the isotropic 
luminosity function, <1>(L) = J d9f b (9)§(L,9), by 
fitting the peak flux distribution using a simple 
parametrization of this function. When consider- 
ing the angular distribution wc also examine (fol- 
lowing Perna, Sari & Frail 2003) the USJ model 
and show that the implied peak flux distribution 
is somewhat inconsistent with the observed one, 
with the problem most severe at the low end of 
the peak luminosity distribution. 

We find that the data support L x 9 2 ~ const 
(see also Van Putten & Regimbau 2003). There 



is a clustering of the true peak luminosity around 
~ 3.2 x 10 49 erg/s, but note that the distribution is 
not a narrow delta-function. We notice then that 
once we have such a relationship the luminosity 
function implies an angular distribution function. 
This is true even when some spread in the value of 
L x 9 2 is taken into account, reflecting the spread 
observed in current data. The resulting angular 
and redshift distributions are consistent with the 
observations. Our analysis is well motivated but 
approximate. Due to width of the luminosity angle 
distribution and given the data quality and selec- 
tion effects, wc believe there is no point in rigorous 
maximal likelihood analysis and the like. We give 
approximate estimates of the model parameters 
and approximate uncertainties. 

The paper is structured in two main parts. In 
the first part we rederive (following Schmidt, 1999) 
the GRB luminosity function. In the second part 
we consider the angular distribution that follows 
from the luminosity function and the peak lumi- 
nosity - angle relation. We compare the predicted 
angular distribution with the observed one. Using 
this distribution we re-calculate the correction to 
the rate of GRBs due to beaming. Unlike Frail et 
al., (2001) we use a weighted average of the angu- 
lar distribution. This yields a significantly smaller 
correction factor. We also estimate a the correc- 
tion factor needed for the USJ model. Note that 
this revised correction factor does not depend on 
the relation between the peak luminosity and the 
angle. 

The observed sample we use to explore the peak 
luminosity - angle relation is small and its selection 
effects arc hard to quantify. Therefore it is difficult 
to give robust conclusions from our analysis at this 
stage. However, future missions and in particular 
SWIFT will allow us, hopefully in the near future, 
to use this procedure and test our conclusions with 
larger and less biased samples. 

2. Luminosity function from the BATSE 
sample 

We consider all the long GRBs (T 90 > 2sec) 
(Kouvcliotou et al., 1993), detected while the 
BATSE onboard trigger (Paciesas et al. 1999) was 
set for 5.5 a over background in at least two detec- 
tors, in the energy range 50-300 keV. Among those 
we took the bursts for which C max / C min > 1 at the 
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1024 ms timescale, where C max is the count rate 
in the second brightest illuminated detector and 
C m in is the minimum detectable rate. Using this 
sample of 595 GRBs we find (U/U max ) = 0.294. 

The method used to derive the luminosity func- 
tion is essentially the same of Schmidt (1999). We 
consider a broken power law with lower and up- 
per limits, 1/Ai and A 2 , respectively. The local 
luminosity function of GRB peak luminosities L, 
defined as the co-moving space density of GRBs 
in the interval log!/ to logi + d\ogL is: 



*o(i) 



(L/L*) a L*/A 1 < L < L* 
(L/L*) 13 L*<L<A 2 L* 



(1) 



where c is a normalization constant so that the 
integral over the luminosity function equals unity. 
We stress that this luminosity function is the 
"isotropic-equivalent" luminosity function. I.e. it 
does not include a correction factor due to the fact 
that GRBs are beamed. 

Following Schmidt (2001) we employ the parametriza- 
tion of Porciani & Madau (2001), in particular, 
their SFR model SF2. 

Rgrb(z) = RsF2{z) 
23exp(3.4z) , ^ ^ . 

where po is the GRB rate at z = and F(z, Om, ^a) = 

[o M (i + zf + n fc (i + zf + r! A ] 1 /2/ ( i + z) v\ 

We consider also the Rowan-Robinson SFR 
(Rowan-Robinson 1999: RR-SFR) that can be fit- 
ted with the expression 

KGRB (z) = po | 10 0.75 Zp eak Z>L ( 3 ) 

For given values of the parameters Ai, A 2 , a, f3 we 
determine L* so that the predicted value (V/V max ) 
equals the observed one and we find the local 
rate of GRB, po, from the observed number of 
GRBs. We then consider the range of parame- 
ters for which there is a reasonable fit to the peak 
flux distribution (sec Figure 1). We use the cos- 
mological parameters H n = 65 km s _1 Mpc -1 , 
Q M = 0.3, and Q A = 0.7. 

The modeling procedure involves the derivation 
of the peak flux P(L,z) of a GRB of peak luminos- 
ity L observed at redshift z: 



P(L,z) 



L C(E 1 {l + z),E 2 (l + z)) 



where Dl(z) is the bolomctric luminosity distance 
and C(E\ ,E 2 ) is the integral of the spectral energy 
distribution between E\ = 50 keV and E 2 = 300 
keV. Schmidt (2001) finds that the median value of 
the spectral photon index in the 50-300keV band 
for the long bursts sample is -1.6 and this can be 
used for a simplified k-correction. We use this 
value for our analysis. 

Objects with luminosity L observed by BATSE 
with a flux limit Pi; m are detectable to a max- 
imum redshift z max (L, Pu m ) that can be derived 
from Eq. 4. The limiting flux has a distribution 
G(-Pii m ) that can be obtained from the distribution 
of C m in of the BATSE catalog. Considering four 
main representative intervals we get that 6% of 
the sample has Pi; m ~ 0.20 ph cm _2 s _1 , 18% has 
P lim - 0.25 ph cm- 2 s-\ 52% has P lim - 0.27 ph 
cm^s- 1 and 24% has ,P lim - 0.32 ph cm^s" 1 . 
The inclusion of this variation of Pu m , and the im- 
plied different samples, is the main difference be- 
tween our analysis of the luminosity function and 
Schmidt's. 

The number of bursts with a peak flux > P is 
given by: 

N(>P) = J^ (L)dlogL J G(P lim )dP lim 
AL - P/P ^ Rgrb(z) dV(z) 



J 

Jo 



l + z 



dz 



dz (5) 



**Di(z) 



C(E\,E 2 ) 



(4) 



where the factor (1 + z)^ 1 accounts for the cos- 
mological time dilation and dV{z)/dz is the co- 
moving volume element. We determine L* so that 
(V/V max ) - 0.294. 

If we approximate Schmidt's results (Schmidt 
2001) the luminosity function can be character- 
ized as two power laws of slopes a = —0.6 and 
= -2, with Ai = 30 and A 2 = 10 and with 
an isotropic-equivalent break peak luminosity of 
L* - 3.2 x 10 51 erg/s and a local GRB p ~ 0.5 
Gpc _3 yr _1 . Using these values we show in Fig- 
ure 1 that the predicted logN-log(P/Pn m ) dis- 
tribution doesn't agree with the observed logN- 
log(C max /C m i n ). In principle one should perform 
a maximum likelihood analysis to obtain the best 
values of the luminosity parameters (and even the 
GRB rate parameters). However, we feel that a 
simpler approach is sufficient for our purpose, es- 
pecially in view of the small size of the sample 
used in the later part of the analysis. We simply 
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vary a and 0, keeping Ai = 30 and A2 = 10 and 
inspect the quality of the fit to the observed logN- 
log(C ma x/C m in) distribution. To obtain the local 
rate of GRBs per unit volume, po we need to es- 
timate the effective full-sky coverage of our GRB 
sample. We find 595 events in this sample that 
were detected over 1386 days in the 50-300 keV 
channel of BATSE, with a sky exposure of 48% . 
We also take into account that this sample is 47% 
of the long GRBs. 

In Figure 1 we show a comparison of the 
observed logN-log(C max /C m i n ) with several pre- 
dicted logN-log(P/Pii m ) distributions obtained 
with the RR-SFR and with a non evolving lu- 
minosity function of the form given in Eq. 1 
(the distributions are similar for the SF2). We 
find reasonable fits with —0.6 < a < —0.1 and 
— 3 < < — 2 taking into account that a and 
are correlated, namely, a from -0.1 to -0.6 require 
decreasing from -2 to -3. As we can see from 
Table 1, among the different acceptable models 
L* varies by a factor of <~ 3. po varies by a factor 
of ~ 2 within a given SFR and as expected it is 
larger by a factor of ~ 2 in the RR model than in 
SF2. 

Since random errors in a cumulative distribu- 
tion like N(> P) propagate in an unknown way, 
we present in Figure 2 our results for the differen- 
tial distributions n(P) = dN/dP. We find reason- 
able fits for the range of a and given above. 

For our analysis we used 595 long GRBs de- 
tected at a time scale of 1024 ms from the BATSE 
catalog. A much larger sample is now available in 
the GUSBAD catalog which lists 2204 GRBs at a 
time scale of 1024 ms 

(http://www.astro.caltech.edu/ mxs/grb/GUSBAD). 
We find that our main results remain the same if 
we consider this catalog. 

From Figures 1 and 2 it is clear that the best set 
of values that reproduce the logN-log(C max /C m i n ) 
distribution is a = —0.1 and = — 2 and these 
will be the values used in the rest of our paper. 
Our results are slightly different from Schmidt's 
mainly because we consider a different sample of 
GRBs and the quantity Pj P; m (which correspond 
to the observed (C max /C m i n )) instead of the flux P 
used by Schmidt (2001). In this way we take into 
account the fact that different bursts are detected 
with different flux thresholds. 



So far we have used a parametric fit for the lu- 
minosity function. Within the USJ jet model, with 
L oc 9~ 2 as forced by the energy-angle relation, 
we have a simple luminosity function: $> S j oc L~ 2 
(Perna et al. 2003). In Figures 1 and 2 we 
also compare the predicted logN-log(P/Pi m ) dis- 
tribution that arises from this luminosity func- 
tion with the RR-SFR with the observed logN- 
log(C max /C m in) distribution. We find that within 
the USJ scenario the distribution is inconsistent 
with the data. The problem arises mostly for low 
peak fluxes. It is important to note, however, that 
our method of GRB selection does not introduce 
a bias against the inclusion of such low peak flux 
bursts. A "quasi-universal" jet (Zhang et al. 2003) 
obviously allows more freedom in fitting the data 
and may be consistent. 

We can use now the luminosity function to de- 
rive the observed rcdshift distribution: 



N(z) 



Rgrb(z) dV{z) 



1 + z 



dz 



A min (P lin 



§ {L)d\ogL, 
(6) 

where I-min(Pim) is the luminosity correspond- 
ing to the minimum peak flux P; m for a burst 
at redshift z and L max = L* x A 2 = 10 L*. 
This minimal peak flux corresponds to the sen- 
sitivity of the gamma-ray burst detector used. 
We use, Pim ~ 0.5 — 1 ph cm _2 s _1 , which is 
roughly the limiting flux for the GRBM on Bep- 
poSAX (Guidorzi PhD thesis). We compare this 
distribution with the observed distribution of all 
the bursts with an available redshift: 32 bursts 
from http://www.mpe.mpg.de/ jcg/grbgcn.html 
(excluding GRB980425 with z = 0.0085). It is 
hard to quantify the selection effects that arise 
in the determination of the redshift. The prob- 
lem is most severe in the range 1.3 < z < 
2.5 where visual spectroscopy (Hogg & Fruchter, 
1999). Following Hogg & Fruchter (1999) we 
consider all the GRB with optical afterglow but 
without a measured redshift to be in this rcd- 
shift range 1.3 < z < 2.5. Using all the bursts 
in http://www.mpe.mpg.de/ jcg/grbgen.html we 
have a sample of 46 GRBs. 

A comparison of the predicted cumulative red- 
shift distribution with the observed one, (see Fig- 
ure 3), reveals that with SF2 we obtain fewer 
bursts at low redshift than the observed ones. This 
happens even after taking into account the selec- 
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SFR 


a 





L* (erg/s) 


Po (GpcrV- 1 ) 


SF2 


-0.1 


-2 


6.3 x 10 bi 


0.18 


SF2 


-0.6 


-3 


1.6 x 10 52 


0.16 


RR-SFR 


-0.1 


-2 


4.4 x 10 51 


0.44 


RR-SFR 


-0.6 


-3 


1.1 x 10 52 


0.19 



Table 1: Limiting values of a and ranges needed to reproduce the observed (C max /C m i n ) distribution and 
the implied values of L* and p for the two SFRs considered in the paper. 



tion effects. In fact a Kolmogorov-Smirnoff test 
(KS) shows that the two distributions are only 
marginally compatible (5% level). The RR-SFR 
seems to reproduce better the redshift distribu- 
tion, as the KS test shows that they are compara- 
ble at the level of 20%. Therefore we will focus on 
this distribution in the rest of the analysis. 

3. Distribution of opening angles 

We turn now to the peak luminosity - angle re- 
lation, which is analogous to the fluence - angle 
relation (Frail et al., 2001; Panaitescu & Kumar 
2001). We consider a sample of 9 GRBs of the 
BATSE 4B catalog for which it has been possible 
to determine the redshift and the jet opening an- 
gle, 9 (Bloom et al. 2003). In order to estimate 
the peak luminosity we need the peak fluxes. The 
advantage in taking only the BATSE GRBs is that 
we have all the peak fluxes estimated in the same 
way (averaged over the 1.024 sec BATSE trigger 
in the energy range 50-300 keV) 

For all the bursts of this sample we estimate the 
peak luminosity along the jet of opening angle 9 
as: 

h = fbL. (7) 

In our analysis we consider the angles up to the 
maximum value implied by the observations (~ 0.7 
rad). Within this range of 9 the beaming factor 
can be approximated as fb = (1 — cos9) ~ 9 2 /2. 
In the upper panel of Figure 4 (left side) we show 
the distribution of Lj obtained from this sample. 
Comparing this distribution with the isotropic one 
shown in the lower panel of the same figure (left 
side) we note that the peak luminosity, when cor- 
rected for beaming, is clustered around a value 
Lj ~ 10 49 5 erg/s. A similar result, shown in Fig- 
ure 4 (right side), is obtained when we consider 
a larger sample, of 19 bursts with an angular es- 
timate from 0.05 to 0.7 radians. Here there are 



larger uncertainties in the determination of the 
peak luminosity, as the bursts were detected by 
different detectors with different averaging time 
and different spectral resolution. We have extrap- 
olated the GRB fluxes to the BATSE range 50- 
300 keV using the method described in Sethi & 
Bhargavi (2001) reelaborated for our spectrum. 
The redshifts and fluxes were taken from the ta- 
ble given in Van Putten & Regimbau (2003) who 
also find that the GRB peak luminosities and the 
beaming factors are correlated. 

We find a correlation coefficient between L and 
1/ fb of ~ 0.5. The probability to obtain such a 
large value by chance, for uncorrelated L and 1/ fb, 
is smaller than 0.05. The regression line slope of 
log(L) vs log(l/f b ) is - 0.8. 

The correlation is supported by a recent discov- 
ery of a relation between the observed peak en- 
ergy, E p , in the spectrum and L (Yonetoku et al. 
2003) . This together with the correlation between 
E p and the isotropic equivalent energy in 7-rays, 
i?i S o, (Amati et al. 2002 and Lamb, Donaghy & 
Graziani 2003) and with the Energy - angle rela- 
tion (Frail et al., 2001; Panaitescu & Kumar 2001) 
imply Lj = fbL. However, within the present sam- 
ple there seem to be some spread in the value of 
Lj as can be seen from Figure 4. To check the 
validity of our analysis we take this spread into 
account. 

We first consider the case in which the distribu- 
tion of Lj is represented by a narrow-delta func- 
tion and then discuss the effect of including the 
observed spread in Lj 

Now if the jet angle and the peak luminosity are 
related we can derive the real distribution of the 
GRB jet opening angles. Let P{9, L) be the dis- 
tribution of number of bursts with opening angle 
9 and isotropic luminosity L. The observed distri- 
bution of number of bursts with opening angle 9 
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in the small angles approximation can be written 
as: 



Pobs(0) 



I 

I 



f b P(9,L)dL~ J 6 lp(9,L)dL 

* ma *(e,p Hm ) Rgrb{z ) dV (z) 



l + z 



dz 



dz (8) 



Given the correlation between the GRB peak lu- 
minosity and the beaming factors we can assume 
that the condition L(9 2 /2) = Lj, with Lj roughly 
constant for all the GRBs, holds. Hence we have 
that P(0,L) = P(L)S(9 - y/2L~/L). Using the 
fact that P(L) — L/Lj$>(L) we obtain the intrin- 
sic (corrected for beaming) angle distribution: 



P{6) 



0*2/30-2,3- 
Q*2a Q—2a- 



9<9* 
9>9* 



(9) 



where c is a normalization constant so that the 
integral over the angular distribution equals unity 
and 9* = \j2LjjL* is the break angle of the 9 
distribution. 9* « 0.12rad with our canonical pa- 
rameters. This equation should be compared with 
P{9) oc sin# ~ 9 that is predicted by the US J 
model. It is clear that the two distributions are 
inconsistent for large angles. A partial agreement 
is possible only for small angles with (3 ~ —2. The 
disagreement at large angles is just another indi- 
cation of the fact that luminosity function implied 
by the USJ jet model is incompatible with the ob- 
served peak flux distribution at the low flux end 
(see Figure 1). 

Substituting this relation in Eq.(8) and after 
some algebra we find the observed 9 distribution: 



Pobs(0) 



Ac L^{2L J /9 2 ) 



I 



9 3 

z max (e,p lim ) 



Rgrb(z) dV{z) 



l + z 



dz 



dz. (10) 



We compare the theoretical distribution with 
a sample of 24 GRBs with measured 9 values 
(Blooom et al., 2003). We stress, again, that selec- 
tion effects are hard to quantify especially for large 
and small opening angles. In Figure 5 we show the 
predicted cumulative 9 distribution derived from 
our luminosity function with a RR-SFR and the 
observed one. We have chosen two values of Lj as 
representative of the distributions given in Figure 
4 (upper panel). We see that with these values of 



Lj the opening angle distribution reproduces quite 
well the observed one. 

We have repeated the analysis with a distribu- 
tion of log 10 (Lj) with finite width, ~ 0.4, reflect- 
ing the spread in Lj values seen in current data 
(Figure 4). The resulting 9 distribution is similar 
to the one that arises from the combination of the 
two values of Lj (black solid line in Figure 5). 

Once we know the angular distribution we can 
turn to our final goal and estimate the true GRB 
rate. To do so we multiply the local observed rate, 
Po obtained in §2 by the correction factor that 
now can be estimated directly from the luminosity 
function given the relation between ft and L: 



(ft 1 ) 



//j- 1 *„(L)d log L fL/Lj$ (L)dlogL 



J$ (L)d\ogL J$ (L)dlogL 

(11) 

ifb 1 ) ^ s the beaming factor averaged over the 
true, intrinsic, distribution of opening angles. We 
expect this average to be smaller that the value 
derived by Frail et al. (2001) due to the follow- 
ing reason: In the intrinsic luminosity distribu- 
tion there are many bursts with a low luminosity 
and large opening angles. These bursts domi- 
nate the rate estimate and this should be factored 
in when calculating the average rate correction 
(Sari, 2003). Put differently, in the observed dis- 
tribution, bursts with large L, i.e. large f b ~ , are 
over-represented compared to the true distribution 
since they are observed to larger distances. In or- 
der to see this explicitly, we note that the (/ b _1 ) 
can be derived directly from the observed distribu- 
tion of GRB peak luminosities without deriving 
$(L), using the following argument. The total 
rate of GRBs (per unit L) in the observable uni- 
verse is dN/d\ogL — po(L / Lj)&(L)Vu , where Vu 
is the effective volume of the observable universe, 
V v = f* m dz{dV/dz)R GRB {z)/R GRB (Q){\ + z) 
(z m is the maximum redshift from which sig- 
nals arrive at us by today). The observed 
rate is dN b s /d\ogL = p $>(L)Vl 7 where Vl 
is the volume out to zl, the maximum red- 
shift out to which bursts can be detected, Vl = 
J a z(L) dz (dV/dz)R GRB {z)/R GRB (0)(l + z). The 
beaming factor, defined in Eq. 11, can therefore 
be written as 



(fb 1 ) 



JdL 



dN obs L V v 
dL L 3 Vl 



f 



dL%*% 



(12) 
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Eq. 12 explicitly reflects the fact that in calculat- 
ing the average beaming factor by which the ob- 
served local rate po should be multiplied in order 
to obtain the true rate, (/r ), the observed dis- 
tribution of bursts should be weighted by Vu/Vl- 
Since Vu/Vl is decreasing with L, a smaller weight 
is given to large L, i.e. to large f^ 1 , than in the 
observed distribution. 

Using Eq. 11 with Lj = 10 49 5 erg/s, we find 
~ 53 for SF2 and (fc 1 ) ~ 60 for RR-SFR. 
These results do not change by inclusion of the dis- 
persion in Lj . As expected, these values are lower 
than the value obtained by Frail et al. (2001), who 
derived (f^ 1 ) by averaging without the oc V^ 1 
weight. Using Eq. 12 with Lj — 10 49 5 erg/s and 
RR-SFR we obtain (J" 1 } ~ 100 for for the BATSE 
sample, and (/ b _1 ) — 50 for the larger sample. Us- 
ing Eq. 12 with Lj determined for each burst from 
its estimated opening angle, we find (f^ 1 ) — 70 
for for the BATSE sample, and (f^ 1 ) ~ 55 for the 
larger sample. The values obtained using these 
different choices are all within (/ 6 _1 ) = 75 ± 25. 
This agreement reflects the fact that our derived 
$(£) gives a fair representation of the data. The 
range of values obtained reflects the range of un- 
certainty given current data. 

We can also estimate the correction factor for 
the rate of GRB in the case of the USJ model 
in the following way: The total flux of GRBs 
per year (or per any other unit of time) is an 
observed quantity that can be obtained by sum- 
ming over the observed distribution. In the USJ 
model the observed energy-angle relation implies: 
E{6) = E o /(tt0 2 ) for 9>9 C and E{6) = E /(ir6 2 c ) 
for 9 < 6 C , where 9 C is the core angle of the jet. 
The total energy that a burst with a USJ emits is: 

£ USJ = 2[ (E /9 2 c ) 2 9d9 + / EaO^M 
Jo Je c 

= £ o [l + 21og(0 max /0c)], (13) 

where 9 max is the maximal angle to which the jet 
extends. 

This immediately implies that 

iVuniform/ATuSJ = [1 + 2 log(9 max /9 c )} (14) 

We don't know for sure what are the upper and 
lower limits but the logarithmic dependancc im- 
plies that the factor cannot be smaller than 2 or 



much larger than 5. This gives us the rate of uni- 
form jets to be about factor of 4 below the number 
of "Uniform" jets. Put differently this suggest that 
the correction factor for the true rate of USJ com- 
pared to the rate of observed GRBs is a factor of 
~ 20 ±10. 

4. Implications to Orphan Afterglows 

The realization that gamma-ray bursts are 
beamed with rather narrow opening angles, while 
the following afterglow can be observed over a 
wider angular range, led to the search for orphan 
afterglows: afterglows that are not associated with 
observed prompt GRB emission. The observations 
of orphan afterglows would allow to estimate the 
opening angles and the true rate of GRBs (Rhoads 
1997). 

Nakar, Piran & Granot (2002) have estimated 
the total number of optical orphan afterglows 
given a limiting magnitude using the true rate of 
GRBs given by Frail et al. (2001). Nakar, Pi- 
ran & Granot (2002) assumed, for simplicity, that 
all bursts have the same opening angle, denoted 
9*. In their canonical model they assume that 
9* = O.lrad while in the "optimistic" model they 
assume 9* — 0.05rad. The smaller opening angle 
gives of course more orphan afterglows. This anal- 
ysis have to be modified now using the new rate 
and the new correction factor that we have found. 
It is clear that the assumption of a fixed opening 
angle is not good enough and to obtain the rate 
of optical orphan afterglows one have to perform 
another weighted average over the observed 9 dis- 
tribution. This weight favors narrow jets which 
produce orphan afterglows over a wide solid angle 
and for which the rate correction is large. Overall 
(Nakar, 2003) one has to correct downwards by a 
factor of 1.6 the rates of the "canonical" model of 
Nakar, Piran & Granot (2002) . The orphan after- 
glow rate obtained in the "optimistic" model are 
clearly an overestimate of the true rate, as this 
model assumed 9* = 0.05 which is significantly 
lower than the average opening angle that we find 
here. 

The situation is different for orphan radio af- 
terglows, which are seen at large angles. We re- 
estimate the number of orphan radio afterglows 
associated to GRBs that should be detected in a 
flux-limited radio survey (Levinson et al. 2002). 
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Levinson et al. have shown that the number of 
such radio afterglows detected over all sky at any 
given time above a threshold /„, m i n at 1 GHz is 

N ~ 10 ^ b ^ f P ° \ ( fv,min \ 

R ~ 70 V0-5Gpc- 3 yr-V VSmJy/ 

3/2 9/8 19/24^11/6 

Here, i? = 10 51 i?5i ergs is the total fireball en- 
ergy, assumed equal for all bursts following Frail 
ct al. (2001) and Panaitescu & Kumar (2001), 
n = 10 _1 n_icm -3 is the number density of the 
ambient medium into which the blast wave ex- 
pands, andes = 10 _1 £b ; _i (e e = 10~°- 5 £ e ,_o.5) is 
the fraction of post shock thermal energy carried 
by magnetic field (electrons). For f^. m - m ~ 5 mJy, 
N R depends mainly on the local rate p and only 
weakly on the redshift evolution of the GRB rate, 
since most of the detectable afterglows lie at low 
redshift z < 0.2 (Levinson et al. 2002). 

Afterglow observations imply a universal value 
of e e close to equipartition, £ e ,-o.5 — 1) based on 
the clustering of explosion energies (Frail et al. 2001) 
and of X-ray afterglow luminosity (Freedman & 
Waxman 2001, Berger et al. 2003) The value 
of eb is less well constrained by observations, 
since in most cases there is a degeneracy be- 
tween eb and n in model predictions that can 
be tested by observations. The peak flux of a 
GRB afterglow seen by an observer lying along 
the jet axis is proportional to E^neB) 1 ^ 2 , and 
for typical luminosity distance of 3 x 10 28 cm 
it is w 10(e B n/10- 3 cm- 3 ) 1 / 2 J B isOi5 4 mJy (Wax- 
man 1997, Gruzinov & Waxman 1999, Wijers & 
Galama 1999). Here, E iso = 10 54 i?i SOi 54 erg is 
the isotropic equivalent energy. Observed after- 
glow fluxes generally imply e^fi > 10~ 3 cm~ 3 for 
-Eiso,54 ~ 0.1, and values esn <~ 10 -1 cm~ 3 are ob- 
tained in several cases. We have therefore chosen 
the normalization n_i = £b,-i = 1 in Eq. (15). 

Using our value for {f b ~ ), (f^~ ) = 70 instead 
°f = 500 given in Frail et al. (2001), reduces 

the number of expected radio afterglows by a fac- 
tor of ~ 7. However, the number of afterglows 
expected to be detected by all sky ~ 1 mJy radio 
surveys is still large, exceeding several tens. 

It should be pointed out that the lower limit on 
the beaming factor inferred from the analysis of 
Levinson et al. (2002) is unaffected by our modi- 
fied value of Assuming a fixed value of fb 



and expressing Nr in terms of the isotropic equiv- 
alent GRB energy, .Ei so . = lO 54 ^^.^ ergs, we 
have 

Nr * W fb U-5G P c- 3 yr-i ) e e,-o. 5 

9/8 19/24^11/6 

Formally, the average value of E^J 6 should appear 
in this equation. We have chosen < E^J 6 > 6 / 11 = 
10 54 erg as a representative value. The upper limit 
on the number of afterglows derived in Levinson 
et al. (2002) implies a lower limit on the beam- 
ing factor, f^ 1 > 40 for the parameter choice in 
Eq. (16), indicating that radio surveys may indeed 
put relevant constraints on the beaming factor. 

5. Conclusion 

In this work a combined L-9 distribution func- 
tion of GRBs is derived, $(!/,#), for the uniform 
jet model. To this aim we have rederived (follow- 
ing Schmidt, 1999) the luminosity function, $(£) 
by fitting the peak flux distribution and used the 
relation L x 2 ~ const implied by the observation 
on the sample considered by Bloom et al. (2003). 
We have compared our results with those obtained 
in the framework of USJ jet model (Perna et al. 
2003) showing that the luminosity function im- 
plied by this model leads to a peak flux distri- 
bution some what inconsistent with the observed 
one. 

The luminosity function that best fits the ob- 
served peak flux distribution is characterized by 
two power laws with slopes a = —0.1 and [3 = — 2 
and isotropic-equivalent break luminosity L* ~ 
7.1 x 10 51 erg/s for a SF2 or L* - 4.4 x 10 51 erg/s 
for a RR-SFR. Repeating the Schmidt analysis we 
have found the observed local rate of long GRBs 
po ~ 0.10 Gpc~ 3 yr -1 and p ~ 0.44 Gpc" 3 yr _1 
for the two SFR respectively. We have also shown 
that with these luminosity functions we find a rea- 
sonable agreement with the observed redshift dis- 
tribution for a RR-SFR, while a SF2 predicts too 
few bursts at low redshift. Nevertheless since the 
redshift sample is strongly affected by selection ef- 
fects that are hard to estimate the possibility that 
GRBs follow the SF2 cannot be ruled out by this 
analysis. 

Using the sample of Bloom et al. (2003) we 
have shown that the energy-angle relation (Frail et 
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al., 2001; Panaitescu & Kumar 2001) applies also 
to the peak luminosity and that there is a cluster- 
ing of the true peak luminosity around ~ 3.2 x 10 49 
crg/s. However the distribution is not a narrow 
delta-function. This implies that a luminosity 
function determines a 6 distribution as L is re- 
lated to 6. This is true even when the spread in 
Lj observed in the current data is taken into ac- 
count. One important result of our work is that 
the resulting angular distribution found under this 
assumption is consistent with the observation. 

We have re-calculated the correction to the rate 
of GRBs due to beaming using a weighted average 
of the predicted angular distribution instead of a 
simple average over the "true" angular distribu- 
tion as done by Frail et al. (2001) and by van 
Putten & Regimbau (2003). We find a correction 
factor 75 ± 25 (see the end of § 3) . This is signifi- 
cantly smaller than the commonly used correction 
factor, ~ 500, estimated by Frail et al. (2001) 
or ~ 475 estimated by van Putten & Regimbau 
(2003). . This correction should also influence the 
estimated rates of both optical (Nakar, Piran & 
Granot, 2002) and radio (Levinson et al., 2002) 
orphan afterglows (see § 4). 
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Fig. 1— The predicted logN-log(P/Pi im ) distri- 
bution for different values of a and (3 with a RR- 
SFR vs. the observed logN-log(C m ax/C m i n ) taken 
from the BATSE catalog. We also plot the pre- 
dicted logN-log(P/Pij m ) distribution in the US J 
jet model, for which the inconsistency at the low 
peak flux range is apparent. 



obs. 

o=-0.1 p=-2 
a=-0.6 p=-2 
USJ model 



I , , , 1 

0.5 1 1.5 2 2.5 

log10(C M /G m ) 



Fig. 2. — The predicted differential distribu- 
tion, n(P/Plim), for different values of a and j3 
with a RR-SFR vs. the observed n(C max /C m i n ) 
taken from the BATSE catalog. We also plot the 
predicted differential distribution in the USJ jet 
model, for which the inconsistency at the low peak 
flux range is apparent. 

This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 3. — The predicted cumulative distributions 
of the GRBs redshift for the two SFR for a = -0.6 
(solid line) and our a = —0.1 (dashed line). The 
histograms show the observed cumulative redshift 
distribution and a distribution where selection ef- 
fects are taken into account assuming that all the 
GRB with no redshift but with optical afterglow 
lie in the range 1.3 < z < 2.5. 
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Fig. 5. — The predicted cumulative distributions 
of GRB opening angles for two different values of 
Lj and for a combination of the two (reflecting 
the spread in the luminosity-angle relation) for 
RR-SFR. These results are compared to the curve 
obtained considering the observed distribution of 
Lj . The histograms show the observed cumulative 
opening angle distribution. 



Fig. 4. — Left: the distribution of the isotropic 
peak luminosity (lower panel) and the distribu- 
tion of the jet peak luminosity (upper panel) of 
10 GRBs in the BATSE catalog with angle deter- 
mination. Right: the distribution of the isotropic 
peak luminosity (lower panel) and the distribution 
of the jet peak luminosity (upper panel) of a large 
sample of 19 GRBs with angle determination. 
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